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Turbulent Taylor-Couette flow between counter-rotating cylinders develops intermittently fluctu- 
ating boundary layers for sufficient counter-rotation. We demonstrate the phenomenon in direct 
numerical simulations for radius ratios rj = 0.5 and 0.71 and propose a theoretical model for the 
critical value in the rotation ratio. Numerical results as well as experiments show that the onset of 
this intermittency coincides with the maximum in torque. The variations in torque correlate with 
the variations in mean Taylor vortex flow which is first enhanced for weak counter-rotation, and 
then reduced as intermittency sets in. To support the model, we compare to numerical results, 
experiments at higher Reynolds numbers, and to Wendt's data. 
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I. INTRODUCTION 

The flow between concentric cylinders has served as a 
paradigm for the transition to turbulence since Taylors 
1923 characterization of the bifurcation from laminar to 
vortical flows Q] . Many transitions between spatially and 
temporally simple flow states like vortices, modulated 
vortices and traveling waves, are accessible by standard 
bifurcation theory and have been described and studied 
in considerable detail [5J [3] . The turbulent states that 
are reached after several of these bifurcations are not al- 
ways homogeneous but can show patchiness in the form 
of turbulent spots or turbulent spirals [3HZ]- In addi- 
tion to such azimuthal and axial modulations, Coughlin 
and Marcus [5] have described a radial inhomogeneity 
for counter-rotating cylinders, which consists of turbu- 
lent bursts that have been seen in experiments by Colovas 
and Andereck [5]. Measurements in the Twente Taylor- 
Couette facility have confirmed the presence of this in- 
homogeneity up to Reynolds numbers of 10 6 [10]. More- 
over, the onset of this instability has been linked to the 
maximum in torque that appears for moderate counter- 
rotation [TOlUT] . 

Coughlin and Marcus [8j already suggested that this 
radial inhomogeneity should be linked to the presence 
of a neutral surface of vanishing angular velocity in the 
laminar Couette profile for counter-rotating cylinders. 
The region between the inner cylinder and the neutral 
surface is inviscidly unstable by the Rayleigh criterion, 
whereas the region between the neutral surface and the 
outer cylinder is Rayleigh stable. These considerations 
do not provide an immediate prediction for the onset of 
inhomogeneity since the neutral surface appears with in- 
finitesimal amounts of counter-rotation already. We here 
discuss the extensions needed in order to derive a predic- 
tive theory for the onset of this intermittency, compare it 
with observations for different radius ratios, and describe 
the link to the torque maxima. 

The outline of the paper is as follows. In section II 
we introduce the numerical simulations and describe the 
phenomenon. We also identify the onset of intermit- 



tency and summarize numerical and experimental results 
for torque maxima, including a reanalysis of the data of 
Wendt [12] which is documented in appendix [Aj In sec- 
tion [HI we present the argument for the boundary layer 



intermittency, and in section IV we describe the link to 



the torque maximum. We conclude with a few remarks 
in section Ivl 



II. BOUNDARY LAYER INTERMITTENCY 

In order to introduce the phenomenon we present re- 
sults from direct numerical simulations for different ra- 
dius ratios. We solve the incompressible Navier-Stokes 
equation with the spectral scheme explained in [13] . For 
the dimensionless units we measure all lengths and times 
in units of the gap width d = r a — r,, where r a and rj 
are the radii of the outer and inner cylinders, respec- 
tively, and the viscous time d 2 /v. To avoid the end- 
effects caused by top and bottom lids in experiments, an 
additional periodicity in axial direction of length L z is 
introduced resulting in an aspect ratio Y = L z /d. Here 
r = 2, which allows for one Taylor vortex pair when the 
outer cylinder is at rest. Fourier modes and Chebyshev 
polynomials are employed for expansions in the two pe- 
riodic and the wall-normal direction, respectively. We 
simulate a domain of reduced azimuthal length, i.e. one 
third for rj — Tijr — 0.5 and one ninth for r\ — 0.71 of 
the full azimuthal length. Consequently the flow field re- 
peats three (nine) times to fill the entire circumference. 
We tested that the shorter azimuthal period does not in- 
fluence the computed torques. The criteria used to test 
and verify the code are detailed in [T3]. Specifically, the 
spatial resolution, characterized by the number of modes 
(N Z) Ntp,N r ) in each direction, is chosen so that three 
convergence criteria are satisfied: torque computed at 
the inner and outer cylinder has to agree within a rela- 
tive deviation of 5 • 10 3 , the expansion coefficients in each 
direction have to cover a range of at least 10 4 and the en- 
ergy dissipation estimated from the torque has to agree 
with the volume energy dissipation rate to within 10~ 2 . 
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All these requirements are met in all simulations shown 
here. 

We consider radius ratios r\ = 0.5 and r\ — 0.71. In or- 
der to assess the influence of the mean system rotation on 
turbulent characteristics, the simulations are performed 
at a fixed shear between the cylinder walls, defined by 
Dubrulle et al. [15] as 

2 

Re s = — — \r/Re - Re^ , (1) 
l + rj 

with Re, = (r Q - r^r^i/v and Re a = (r D - ri)r w /v. 
Here, Wj and uj denote the angular velocity of the in- 
ner and outer cylinder and v the kinematic viscosity. For 
both radius ratios we realize various mean rotations char- 
acterized by the rotation ratio fi = uj /uii for the same 
shear Res = 2.0 • 10 4 . According to Lathrop et al. [TS] 
this value is high enough that the vortices noticeable for 
lower Re have disappeared for fj, = 0. The dimension- 
less torque G exerted on the inner and outer cylinder is 
obtained from 

G(t) = v~ 2 r = v-\* ((u r u>) A{r) - vd r (u) A(r) ) , (2) 

where u r and uj = u v /r denote the radial and angular 
velocity, respectively, and (• • -) A r r \ stands for an area av- 
erage over the surface of a concentric cylinder |17) . The 
long-time mean value is obtained from an additional av- 
erage over time. 

For co-rotation, i.e. /i > , torque values at the inner 
and outer cylinders agree on average and in their fluctua- 
tions [13] • However, the situation drastically changes for 
strong counter-rotation as illustrated by the torque time- 
series for /i = —0.5 in Fig. [IJa). Near the inner cylinder, 
fluctuations are small and without long time variation. 
At the outer cylinder the torque exhibits fluctuations of 
relatively strong amplitude and slow dynamics that qual- 
itatively differ from the ones at the inner cylinder. These 
slow fluctuations reflect an intermittent turbulent activ- 
ity in the vicinity of the outer cylinder as demonstrated 
by the cross-flow energy in Fig. 1(b). The cross-flow 
energy, 

E cf (r,t) = (ul + ul) A{r) (3) 

measures the energy content in the transverse velocity 
components at a radial distance r and at an instant of 
time. Small values indicate a flow close to laminar, large 
values a turbulent state. One notes that near the inner 
cylinder the flow is more or less homogeneously turbu- 
lent, whereas towards the outer cylinder one observes an 
increased activity synchronized with the increase in mean 
torque. 

For this geometry, the Rayleigh-stable region 
of the laminar profile extends over the interval 
(r — fj) / d € [0.414, 1], with the lower bound marked by 
a dashed line in Fig. [ljb). The presence of the bursts 
shows that this stability is not maintained, but the 
position of the laminar line still marks the transition 
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FIG. 1: Turbulent bursts for 77 = 0.5, Re s = 20000 and 
li = —0.5. (a) Time-series of torque at the inner cylin- 
der (blue) and the outer cylinder (red) divided by its lam- 
inar value, (b) Spatio-temporal plot of cross-flow energy, 
u 2 + u 2 z \ , averaged along the periodic coordinates with 

/ A(r) 

red (blue) indicating high (low) energy. The dashed line in- 
dicates the laminar neutral surface. 

between small radial variations near the inner cylinder 
and larger ones towards the outer cylinder. 

For Res ~ 2.3 • 10 3 Coughlin and Marcus [8] described 
the bursts as instability of spatially ordered 'interpene- 
trating spiral' flow. We here observe turbulent flow in 
the Rayleigh unstable inner region at Res — 2.0 • 10 4 
that is accompanied by intermittent turbulence in the 
outer one. The connections to spirals cannot be followed 
up because of the limited radial and azimuthal domain 
size that can be computed. The experiments of [10] ex- 
tend this observation to even larger Reynolds numbers 
Re s ~ 10 6 where the presence of turbulent bursts in the 
outer region is deduced from a bimodal distributions of 
the angular velocity. 

The onset of the intermittent behavior is accompanied 
by an increase in the torque fluctuations. Therefore, we 
study the standard deviation oq of the torque relative to 
the mean for the outer cylinder, i.e. the ratio ac/G, as an 
indicator and deduce a critical value fJ, c {f]) for the onset 
from the requirement that aa/G exceeds the base level 
for (j, = 0, cf. Fig. [2] The choice of a G /G as measure 
for the transition is supported by the observation that 
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FIG. 2: Standard deviation oq of temporal torque fluctua- 
tions (divided by the mean G) computed at the inner cylin- 
der (blue circles) and the outer cylinder (red triangles) for a 
constant shear and various global rotations. The dotted line 
marks the fluctuation base level, the dashed line a linear fit 
to the data and the solid line the intersection point /j, c . (a) 
rj = 0.5 and Re s = 20000. (b) 77 = 0.71 and Re s = 19737. 
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FIG. 3: Dependence of the total torque on the global rotation 
for 77 = 0.5 and Re s = 20000 (blue circles) and for r? = 0.71 
and Res = 19737 (red triangles). Error bars based on the 
statistical uncertainty would be smaller than the symbols. 



ing Res [H] • Figure [3] shows the computed torques for 
Res = 2.0 • 10 4 just at the beginning of the asymp- 
totic regime. We determine the rotation ratio of optimal 
transport as the maximum of a quadratic fit, G/G\ am — 
C2fJ? + ci/i + cq, to five data points and find 



Mmax 
A^max 



(0.5) 
(0.71) 



-0.195 ±0.019 
-0.357 ± 0.060. 



(5) 



The uncertainties are deduced from the relative confi- 
dence interval AG/G which results from temporal torque 
fluctuations. This uncertainty in the torque values trans- 
forms into an uncertainty in the maximum location of the 
quadratic fit, i.e. 



it remains relatively unaffected by variations of 11 at the 
inner cylinder. For the numerical simulations we find the 
critical values 

/i c (0.5) = -0.207 ±0.014 
/i c (0.71) = -0.351 ±0.050, (4) 

marked by vertical lines in Fig. [2] The uncertainties 
are estimated as half the gap in fi between the two data 
points next to the maximum. The critical point varies 
with 77, and we will return to this feature in section |III| 
The onset of intermittent bursts was determined experi- 
mentally at /i c (0.716) « -0.368 for Re s ~ 10 6 10,. The 
differences between this value and our observation Q 
may be due to our uncertainty in /_i c or the difference in 
Reynolds number, but they are not significant. 

The rotation ratio /J m ax('7) of maximal torque was 
identified independently in two experiments, (TO] IHJ 
118] . For a constant shear, they report a torque max- 
imization for /j max (0.716) = -0.33 ± 0.04 [TU] and for 
Mmax(0.7245) = —0.333 [11] . Numerical simulations show 
that the torque maximum for counter-rotation only oc- 
curs after a shift of the maximizing /1- value with increas- 



A/i max - J --77--Q ( G / G lam) max (6) 

with the fit coefficient ci and the maximal rescaled torque 

\G J (jlam)max- 

For both radius ratios, the maximizing global rotation 
/Anax compares well with the transition to the radial in- 
homogeneity at fi c , as it was previously observed exper- 
imentally by van Gils et al. [TUJ for 77 = 0.716. On the 
other hand, the critical values for rj = 0.71 are more un- 
certain, which complicates the identification of a corre- 
spondence between them. We note that the torque max- 
imization at /i ma x(0.71) = —0.357 falls in line with the 
experimental observations /j max (0.716) = —0.33 ± 0.04 
and /j max (0.7245) = -0.333 (23111]. 

A final data point is provided by Wendt's data [12], 
which are reanalyzed in the way used in the recent ex- 
periments in appendix [X] Wendt worked with the radius 
ratio of r\ = 0.680, and the maximum in his torque data 
lies near 

Mmax(0.680) = -0.295 ± 0.113 . (7) 
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III. ONSET OF THE RADIAL 
INHOMOGENEITY 

In order to connect the onset of fluctuations with the 
rotation ratio for maximal torque and its dependence on 
the radius ratio, we ask the following questions: (i) What 
is the physical mechanism that determines the onset of 
fluctuations and can one derive a prediction for the onset 
from it? (ii) Why should the rotation ratio of torque 
maximization coincide with the onset of fluctuations? 

A first answer to these questions is provided by van 
Gils et al. [TD]: They suggest that the rotation ratio 
for maximal torque is determined by locations in param- 
eter space (Rei, Re a ) that are equally distant from the 
Rayleigh stability lines fi — r] 2 and /i = — oo. This con- 
dition results in the so called "angle bisector" , 



V 



He 



/^max /-^pred l^bis 



MLSC 



-n 



tan [f — | arctan(ry" 



(8) 



for the location of the torque maximum |10) . Sec- 
ondly, they argue that the onset of fluctuations has 
to coincide with ^ since the intermittent behavior 
in the outer layer reduces radial transport of mo- 
mentum the torque starts to drop. While the an- 
gle bisector agrees well with their measured optimum 
A t max(0-716) = —0.33, it disagrees with our simulation re- 
sult for r\ = 0.5, since /J, ma ,x{0-5) — —0.195 whereas Eq. 
(§) gives /Lt bls (0.5) = -0.309. 

We here propose an explanation for the onset of fluc- 
tuations that is not based on the stability of laminar 
flow but on properties of turbulent flows in general and 
turbulent Taylor vortices in particular. The key idea is 
that the turbulent flow detaches because the inner part 
that is driven by the Rayleigh-unstable region is not suf- 
ficiently strong to maintain persistent turbulence across 
the Rayleigh-stable region to the outer cylinder. But 
the outer region cannot return to laminar for all times, 
because the turbulent transport and the friction has to 
be the same, independent of the radial distribution [T7] , 
The radial range over which the inner unstable region can 
maintain a turbulence is somewhat larger than the inner 
Rayleigh-unstable region. The simple inviscid stability 
calculations for counter-rotation alluded to before gives 
a neutral surface at radius 



1 



T] A — fJL 



(9) 



that separates stable from unstable flow [19] and im- 
plies a detachment of the unstable flow for any fi < 0. 
However, experiments and viscous calculations show that 
Taylor vortices extent beyond this neutral surface when 
counter-rotation sets in, see for example pQ. Esser and 
Grossmann |20j deduced from their stability calculation 
that flow structures protrude the neutral surface by a 
factor 0(77), i.e. the effective extension of secondary flow 
is 



0.5 -0.207 
0.68 

0.71 -0.351 
0.716 -0.368 
0.7245 



-0.195 
-0.295 
-0.357 
-0.33 
-0.333 



-0.191 
-0.321 
-0.344 
-0.349 
-0.356 



-0.309 -0.223 
-0.360 

-0.367 -0.357 

-0.368 

-0.370 



TABLE I: Rotation ratio of the onset of intermittency \i c 
and of the torque maximum /J max together with the new pre- 
diction /ip rc d ( 12 1 and the angle bisector [ibis Q for various 
radius ratios. The rotation ratio /jlsc of the maximal mean- 
flow contribution to the torque is also given for the numerical 
simulations in the last column. 



with a factor 



i( V ) = (1 - 77) 



1 (1 + 77) 3 

2(1 + 3rj) 



(11) 



that takes values between 1.4 and 1.6. The rota- 
tion ratio /x pre( j where the unstable flow detaches from 
the outer cylinder wall then follows from the condition 
r eg (Mprcd) — r Q and reads 



Mpredfa) = -IT 



(a 2 - 2a + 1)77 



- 1 



(2a -1)77 + 1 



(12) 



rEc(p) = n + 0(77) (r„ - n) 



(10) 



For /i < /jp r cd the turbulence can no longer fill the whole 
cylinder gap, i.e. teg < r o, and intermittency has to 
set in. The extended range of the inner unstable region 
as captured by the factor 0(77) then mandates a mini- 
mal counter-rotation for the neutral surface to fall inside 
the cylinders. The evaluation of (12 1 yields predictions 
AVod(0.5) = -0.191 and /Jprcd(0.7l7^= -0.344 that com- 
pare well with the empirically found onsets of intermit- 
tency Q, see also Tab. [I] 



IV. ENHANCED LARGE SCALE 
CIRCULATION 

We now turn to the question of why the torque maxi- 
mum coincides with the onset of the radial intermittency. 
Van Gils et al. c 10| argue that the torque decreases when 
turbulent bursting sets in because of the reduced radial 
transport. However, to obtain a torque maximum at 
the bursting onset fi c , additionally, the torque has to 
increase with /i decreasing from zero to fi c < 0. We 
argue that this increase is caused by a strengthening of 
the mean Taylor vortex flow. Such a large-scale circu- 
lation (LSC) is able to effectively transport momentum 
and, thus, to increase the torque in addition to turbu- 
lent fluctuations. The strengthening of the LSC is due 
to a change in the effective outer boundary condition: 
While for fi > /i prc d the LSC seeks to extend beyond 
the outer cylinder and is restricted by the rigid wall, i.e. 
teg > r o, the rigid boundary conditions become replaced 
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contribution. Moreover to accurately capture the mean 
Taylor vortex motion, we axially shift the instantaneous 
flow fields during the temporal average of u and in ( 14 ) 
so that Taylor vortices always stay at a fixed height. 

For \i > 0, the torque is mainly caused by turbulent 
fluctuations and the mean-flow contribution nearly drops 
to the laminar level, as shown in Fig. [4] Turbulent 
fluctuations also dominate the torque for strong counter- 
rotation, i.e. /i — 0.5 for r\ = 0.5 and fi < 0.71 for 
j] = 0.71. For intermediate rotation ratios, mean-flow 
vortices (LSC) contribute the major share to the torque. 
Note that the onset of mean vortical flow for fi > —r), 
with fi = —r\ corresponding to perfect counter-rotation 
Re Q = — Rei, was previously observed by Ravelet et al. 
[21] . Additionally, the LSC contribution grows with /x de- 
creasing from zero which is consistent with our picture of 
a change in the outer boundary condition from no-slip to 
a less restrictive free-surface condition. The mean Taylor 
vortices are strongest (as measured by their contribution 
to the torque) at the rotation ratio hlsc(v) where G is 
maximized. Using a quadratic fit with an uncertainty 
estimation analogue to ^ we find 



Mlsc(0.5) 
Mlsc(0.71) 



-0.223 ±0.018 
-0.357 ±0.075. 



(15) 



FIG. 4: Rotation dependence of the mean-flow (blue circles) 
and turbulent contribution (red triangles) to the torque (open 
circles). The torque due to the mean flow exhibits a maximum 
marked by the vertical gray line. Dotted lines connecting the 
data serve as guide to the eye. (a) r\ = 0.5 and Res = 20000. 
(b) V = 0.71 and Re s = 19737. 



by a softer free-surface-like laminar outer boundary layer 
as fi approaches /i prc d ~ Mc- Less restricted, the large- 
scale vortices can become stronger before they will be 
destroyed by the bursting for /x < fi c . 

To test this picture we decompose the flow into the 
LSC contribution u = (u) t that contains the mean 
variations in radial and axial direction and the turbu- 
lent fluctuations u' = u — u. With this decomposition 
one finds from Eq. (§ a partitioning of the torque 



G = G + G' (13) 
with the mean-flow (LSC) and turbulent contribution 

G = »~ 2 (r 3 (^u r Q} A(r)!t -vd r {Q) A(r) ^) r , 
G' = v-i(r 3 (uW} A(rht ) r . (14) 

The mixed terms (u r u') and (u' r Q) in ([2]) as well as 
{oj') a ,s t vanish due to the definition of u. While the 
complete torque G is radially independent, similar ex- 
pressions for mean-flow and turbulent contribution vary 
with the radius. Therefore, we introduced an additional 



radial average in ( 14 1 to measure the mean weight of each 



Consequently, the rotation ratio of optimal momentum 
transport by the mean flow coincides with the empirically 
found onset of intermittency, Eq. Q, within the given 
uncertainties. Furthermore, the mean-flow contribution 
is responsible for the maximum in the total torques, cf. 
Fig. [1] and Eq. ([5| , thereby establishing the connection 
between the onset of intermittency and the torque max- 
imum within the framework depicted above. 

This connection implies that the prediction for the in- 
termittency onset also acts as a prediction for torque 
maxima. We thus compare the predictions from the 
boundary layer argument /K pr ed(^) from Eq. (12) and 
the bisection argument /ibis (jf) from Eq. (8) |10j with ex- 
perimental and numerical results for torque maxima in 
Fig. [5] (and Tab. [[]). The rotation ratio /Lt = -2.797 of 
the turbulent bursts found by Coughlin and Marcus [8] 
lies below both predictions and, thus, clearly in the in- 
termittent range. Our simulation result for r\ = 0.71 and 
Wendt's experimental result for r\ = 0.680 are consistent 
with both the angle bisector and the current prediction 
within the error bars. Moreover, we note that the torque 
maximum /x max (0.7245) = —0.333 measured by Paoletti 
and Lathrop [TT] as well as /i majc (0.716) = -0.33 ± 0.04 
measured by van Gils et al. [10] tend towards our pre- 
diction, Eq. p^j ). However considering the usual error 
bars, these values are also consistent with the angle bi- 
sector line. The simulation results for r\ = 0.5 provide a 
better test of both predictions, since the values obtained 
from (8) and (12) differ. The numerical data are in better 
agreement with the boundary layer estimate (12). 
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FIG. 5: Location of torque maxima for different radius ratios 
77. Blue: Angle bisector (Jsj> . Red: boundary layer estimate 
(12 1. The black circle, purple triangle and green star show 
torque maxima from the experiments |1UL|11] and [T2] . The 
cyan squares mark the simulation result J5J. 



V. FINAL REMARKS 

The analysis presented here support the idea that the 
torque increases with increasing counter-rotation because 
the vortices gain in strength until they can no longer sus- 
tain turbulence all across the gap. The torque drops for 
stronger counter-rotation, when the detachment of mean 
vortices from the outer layer leads to radial intermittency. 

The predictions from the boundary layer argument 
presented here and the angle bisection proposal of van 
Gils et al. [TU] give indistinguishable predictions for 
a radius ratio of 77 w 0.75, but the shape of the 77- 
dependencies is sufficiently different that data in particu- 
lar for smaller 77 should allow to distinguish between the 
two. The available data for 77 = 0.5 are in good agreement 
with the present argument. Clearly, results for more 77 
are required and work along those lines is in progress. 

For larger 77 in the limit 77 — > 1 the current theory pre- 
dicts a maximum for a rotation ratio of —0.6 , whereas 
the angle bisection gives a value close to —0.4. How- 
ever, this limit is delicate because the linear instability 
disappears and both theories will most likely have to be 
refined or replaced. Evidence for this is provided for in- 
stance by the measurements by Ravelet et al. [ST] for 
77 = 0.917, which do not show a torque maximum for 
counter-rotation. Figure [5] reflects this uncertainty in the 
prediction by the change from a continuous to a dashed 
line for 77 > 0.9. 

Furthermore, numerical simulations revealed that at 
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FIG. 6: Torques measured by Wendt [12] for 77 = 0.680 and 
various rotation ratios. The values are compensated with the 
laminar Gi am and with the effective scaling G/Giam ~ Res ' 7 
reported by Wendt. The torques for the range shaded in gray 
are further analyzed in Fig. [7] . 



lower Res < 4 • 10 3 the torque is maximized at /i « for 
77 = 0.71 [2]. This larger rotation ratio is not covered 
with the current theory, so that further refinements are 
needed for lower Reynolds numbers and mildly turbulent 
flows. 
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Appendix A: Re-analysis of Wendts data 

Recent experimental studies analyze the dependence 
of torque on the shear rate and on the mean system 
rotation independently. This decomposition is advanta- 
geous since torques can be compensated either by di- 
viding by the effective scaling with the shear [18] or by 
talking the ratio to G(p = 0) llj to study the rota- 
tion dependence. The resulting torque amplitudes are 
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FIG. 7: Compensated torques by Wendt for r\ — 0.680 aver- 
aged in the range 7.8 ■ 10 4 < Res < 1.3 • 10 5 for each rota- 
tion ratio independently. The blue line indicates a quadratic 
least-square fit to the four largest values. Its maximum 
Mmax = —0.295 is marked by the dashed line. 



from figure 9 in [12]. Figure [6] shows the torques for 
various rotation ratios compensated with Res°' 7 which 
Wendt found as effective scaling for 10 4 < Re s < 10 5 . 
One easily sees that the torque depends on the mean ro- 
tation with the largest values for high Res at fi = —0.25. 
We closely follow the analysis in [lOl ITT1 ITS] and average 
the compensated torques in the range 7.8 • 10 4 < Re s < 
1.3-10 5 to find amplitudes depending on the rotation only, 
see Fig. [7] We chose a Reynolds number range that starts 
after the shift of the torque maximum |14j and includes 
the highest data points for —0.50 < fi < —0.17 (cf. Fig. 
[6|. One observes a maximum in the statistical more sig- 
nificant mean amplitudes for moderate counter-rotation 
which was also found in recent studies [TQl HU HE] and 
in current simulations. Based on a quadratic fit to the 
largest amplitudes we find 



Mmax(0.680) = -0.295 ± 0.113 . 



(Al) 



based on numerous measurements for each rotation ra- 
tio which improves statistical significance. In contrast, 
Wendt presented the torque-dependence on the rotation 
for some selected shear Reynolds numbers in figure 10 of 
Ref. |12j . Since this evaluation is based on single mea- 
surements, uncertainties may play a major role. 

Therefore, we here apply the current analysis method 
to Wendt's torque measurements for r/ = 0.680 digitized 



with the uncertainty calculated in analogy to Eq. (J6J) . 
Its relative high level is due to the broad maximum 
in Fig. [7] and due to the few rotation ratios investi- 
gated by Wendt. In spite of the high uncertainty, the 
torque maximization for counter-rotation, i.e. fi max < 0, 
is clear without ambiguity. Moreover, the new maxi- 
mum, /i max (0.680) = —0.295, lies consistently between 
the maxima identified here, cf. Eq. |5|. 
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